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ABSTRACT 

We present a photometric and spectroscopic study of the white dwarf population of the populous, 
intermediate-age open cluster M35 (NGC 2168); this study expands upon our previous study of the 
white dwarfs in this cluster. We spectroscopically confirm 14 white dwarfs in the field of the cluster: 
12 DAs, 1 hot DQ, and 1 DB star. For each DA, we determine the white dwarf mass and cooling age, 
from which we derive the each star's progenitor mass. These data are then added to the empirical 
initial- final mass relation (IFMR), where the M35 WDs contribute significantly to the high-mass end 
of the relation. The resulting points are consistent with previously-published linear fits to the IFMR, 
modulo moderate systematics introduced by the uncertainty in the star cluster age. Based on this 
cluster alone, the observational lower limit on the maximum mass of white dwarf progenitors is found to 
be ~ S.IMq — 5.2Mq at the 95% confidence level; including data from other young open clusters raises 
this limit as high as T.IM©, depending on the cluster membership of three massive WDs and the core- 
composition of the most massive WDs. We find that the apparent distance modulus and extinction 
derived solely from the cluster white dwarfs {{m-M)v = 10.45±0.08 and E{B-V) = 0.185 ± 0.010, 
respectively) is fully consistent with that derived from main-sequence fitting techniques. Four M35 
WDs may be massive enough to have oxygen-neon cores; the assumed core composition does not 
significantly affect the empirical IFMR. Finally, the two non-DA WDs in M35 are photometrically 
consistent with cluster membership; further analysis is required to determine their memberships. 
Subject headings: white dwarfs — open clusters and associations: individual (M35) — stars: evolution 
— supernovae: general 



1. INTRODUCTION 

Intermediate-mass stars {M ~ 6Mq — IOA^q) end their 
evolution in one of two ways, either exploding as a core- 
collapse supernova or losing large amounts of material 
to form a massive white dwarf (WD) star. Recent mod- 
els of massive asymptotic giant branch (AGB) stars sug- 
gest that, for stars with initial (i.e., zero-age main se- 
quence) masses < 7.25M0, the endpoint of stellar evo- 
lution will be a carbon-oxygen (C/0) WD. For stars 
with initial masses ^ 7.25Mq — 9.0Mq, carbon will burn 
in a "super-AGB" star, forming a degenerate oxygen- 
neon (ONe) core which could become an ONe WD (e.g . , 
iGarcia-Berro et all Il997t iPoelarends et al.ll2006L 120081 ). 
More massive super-AGB stars may explode as electron 
capture supernovae, and stars with masses > IIMq ex- 
plode as canonical core-collapse supernovae. 
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The initial mass dividing supernova progenitors from 
white dwarf progenitors, Mup, (also known as and 
referred to as Merit in our earher work) is therefore likely 
to lie in this mass range (~ 7Mq — 9Mq); indeed, cur- 
rent super-AGB models predict M^p ^ 8 — 9Mq, de pend- 
ing o n metallicity and the degree of overshooting (jSiessI 
|2007( ). However, observational constraints on Mup have 
been slow in comin g. The best published observations 
give Mup = 8+2Mg (jKoester fc Rcimers 1996). 

Tight constraints on Mup are important; for a burst 
of sta r formation wit h a Salpeter initial mass function 
(IMF; lSalpeteilll955| ). the number of stars in the mass 
range of 6M0 — IOMq is equal to the number of stars 
with M > IOMq. This factor of two uncertainty in the 
number of supernovae in a starburst region has impor- 
tant implications for quantifyi ng supernova feedback o n 
the interstellar medium (e.g., iMcKee fc Ostrike^ll977f) . 
and for understanding supernova-driven winds in galax- 
ies (e.g.,lbekel & Silk 1986; Martin 2005). 

Observations of WDs provide a lower limit on the value 
of Afup- In the simplest form, this limit can be de- 
termined by identifying WDs in progressively younger 
star clusters until no WDs are found, indicating that the 
turnoff-mass stars are going supernova instead of form- 
ing WD remnants. Variations on this method were used 



2 



Williams, Bolte, & Koester 



by iRomanishin &: Angell (|1980f ) and I Anthony- Twarogj 
(|1982f ). who independently determined Mup > 4Mo — 
5Mq. This method is improyed upon by using spec- 
troscopic analysis of the WDs in each cluster to deter- 
mine each WD's progenitor mass. The most massiye WD 
progenitor is then a lower limit on M,,p. S uch analysis 
was performed by lKoester fc ReimersI ( 1996( ) in the open 
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cluster NGC 2516, where they found M, 
In our initial study of WDs in M35, we determined 
Mup > 5.8Mq (iWilhams et al.li2004t ). 

The relationship of a WD's mass to that of its pro- 
genitor, the initial-final mass relation (IFMR), is also a 
matter of keen interest. This relationship quantifies the 
integrated mass lost by a star oyer its entire eyolution, 
and is therefore a necessary part of understanding chem- 
ical e nrichment and star f ormation efficiency in galaxies 
(e.g., iFerrario et al.l [2005f l. The IFMR represents one 
of the best obseryational constraints on AGB star mass 
loss. These stars are thought to be the primary sites 
of s-process production and to play crucial roles in the 
abundance ratios of helium, ca r bon, arid nitr ogen (e.g., 
iBusso et al.lll999l:lMarigoll200lUHerwidl2005l and refer- 
ences therein). 

Due in large part to the steepness of the IMF, the 
IFMR has relatiyely few points in the high initial mass 
(Afinit > 4M0) region. Four massi ye WDs are known 
from the open cluster NGC 2516 (jKoester fc ReimersI 
Il996l ). The Pleiades proyides one to three points, 
depending on whether the massiye WDs GD 50 and 
PG 01 36-1-251 were once cluster members (jDobbie et al.l 
l2006bD . Age-dating of Sirius A has allo wed the pro- 
genito r mass of Sirius B to be calculated (|Liebert et al.l 
l2005bD . A new paper by Dobbie2008 identifies three mas- 
siye WD progenitors in NGC 3532 and two in NGC 2287. 
Finally, the older open cluster NGC 2099 has one high- 
mass WD, albeit with yery large error bar s on both the 
initial and final mass (iKalirai et al.ll2005b[ ). 

It is suspected that the IFMR should haye some metal- 
hcity dependence (e.g., iMarigo fc Girardil 120071 ). with 
more metal-rich systems producing lower mass WDs. 
Indeed, the WDs in the super-met al-rich cluster NGC 
6791 haye surprisingly low masses (jKalirai et al.ll2007h . 
though the inyoked enhanced mass loss due to the 
high metallicity results in stars circumyenting the AGB 
phase altogether. Comparison of the WD masses in the 
open star clusters NGC 2099, the Hyades, and Prae- 
sepe may show some eyidence of metallicity dependence 
(jKalirai et al.l[20"05b[ ). though this claim hinges on a yet- 
unpublished significantly sub-solar metallicity measure- 
ment for NGC 2099. Published spectroscopic mctallic- 
ity measureme nts for this cluster (Marshall et al. 2005; 
iHartman et al.l 2008) claim a solar metallicity or slightly 
super-solar metallicity for NGC 2099; the younger clus- 
ter age resulting from a higher metallicity would erase 
most of the claimed signal. In short, the metallicity de- 
pendence of the IFMR, while expected, has yet to be 
obseryed conclusiyely. 

Addressing the questions of M^p, the shape of the high- 
mass end of the IFMR, and any metallicity dependence 
yia obseryations thus requires WD obseryations in mul- 
tiple star clusters that are young (ages '--^ 50 — 250 Myr), 
relatiyely rich, and of markedly different metallicities. 



1.1. The Open Cluster M35 

The open cluster NGC 2168 (M35) is an ideal labora- 
tory for addressin g these issues. The cluster has an age o f 
~ 150-200 Myr (jSung fc Bessel]lfT99irwm HiDDe]|[2005h 
and a significantly sub-solar meta llicity [Fe/H]Ri —0.2 
(|Barrado y Nayascues et aLll2001a[ ): newer spectroscopic 
abundance measurements from the WIYN Open Cluster 
Study confirm the sub-solar metallicity ([Fe/Hjw —0.14; 
A. Steinhauer, priyate communication). For the rest of 
this paper, we adopt [Fe/II]= —0.2. The cluster main- 
sequence turnoff mass for the cluster is 4.OM0 — 4.6M0 
based o n the most recent Padoya stellar evolutio nary 
models (|Marigo fc Girardill2007l : iMarigo et all l2008V All 
cluster WDs therefore are the remnants of some of the 
most massiye stars that form WDs. 

Several determinations of the cluster distance mod- 
ulus and extinction exist in the literature. For the 
sake of comparison, we have adjusted these published 
va lues to our assumed met allicity using the relation 
of iPinsonneault et al.l (|1998f ). Resulting apparent dis- 
tance moduli (m — M)v include 10.3 ([Twarog et al.l 
I199I. 10.4 ± 0.1 dSrn^ fc Besselll [T 999l). 10.42 ± 0.13 
(IKalirai fc Tosill2004D . and 10.26 ±0.12 (S araiedini et all 
I2004D . Published values of the r eddening span color ex- 
cess values of E(B -V) = 0.19 (iTwarog et al.l [l99l to 
E{B-V) = 0.255 (iSung fc Besselllll999D . For this pa- 
per, we adopt (m - M)v = 10.3 ± 0.1 and E{B-V) = 
0.22 ± 0.03 as representative of the magnitude of and 
scatter in these values. 

The first search for WDs in M35 was undertaken 
by IRomanishin fc Angel! (|1980 l). who used photographic 
images of the cluster to identify four WD candidates. 
iReimers fc Koesterl (jl98 8') obtained spectra of three of 
these stars, confirming that all three are WDs and es- 
tablishing two as likely m embers. In recent years, deep 
CCD imagin g stud ies of Ivon Hippel et ai] (|2002l ) and 
IKalirai et al.l (j2003D identified several faint, blue stars 
likely to be WDs, but these studies lack the spectroscopic 
observations necessary to determine the WD masses and 
cooling ages. 

Therefore, we began our own imaging and spec- 
troscopic study of t he WD population of M35. In 
[Williams et al.l (j2004f ). we presented initial data on eight 
DA (hydrogen-atmosphere) WDs in the field of NGC 
2168, seven of which were claimed to be cluster mem- 
bers. In this paper, we present data for an additional six 
WDs, four of which are DAs, as well as significant addi- 
tional data on three of the previously-published DAs. All 
of the older data have been completely re-reduced to cor- 
rect some errors in the initial reduction, to incorporate 
improvements to the atmospheric fitting routines, and to 
use the newest WD and stellar evolutionary models. The 
fits in this pap e r ther efore supersede those published in 
iWilliams et all (I200I . 

2. CLUSTER PHOTOMETRY 

2.1. Observations 

Initial UBV imaging of M35 was obtained on 2001 
September 22 UT using the Prime Focus Camera (PF- 
Cam) on the Lick Observatory 3m Shane telescope. 
Weather was photometric, and seeing was moderate and 
steady at « 1'.'6. The imaging covered two ~ 10' x 10' 
fields, one centered ~ 4f5 northeast of the cluster center. 
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and the other ~ 4.5 southwest of the cluster center. As 
data from earher in the ni ght showed evideiice for a non- 
Hnear CCD response (see iwilhams fc Boltell2007l ). these 
photometric data are not presented in this paper. How- 
ever, they were of sufficient quaUty to identify blue-excess 
objects as WD candidates for initial spectroscopy. 

Additional UBV imaging of the cluster was obtained 
on 2004 January 23 UT with the Mosaic- 1 camera on the 
Kitt Peak Mayall 4m telescope. Weather was photomet- 
ric and seeing was excellent: 0'.'75 in V, O'.'Q in B, and 
1'.'4 in U. The Mosaic camera, which has a field-of-view 
of 36' X 36', was centered on the cluster core. In each 
filter, three long exposures were obtained for a total of 
2700s in U and 360s in B and V. Between each exposure, 
a dither of ^ 52" north and east was executed. Single 
shorter exposures of 30s and 5s were also obtained in 
each filter. 

2.2. Photometric reduction 

The Mosaic data were reduced using the MSCRED 
package of IRAF.'^ We closely followed the prescription 
formulated for the NOAO Deep Wide-Field Surv ey (ND- 
WFS) described in online notes bv lJannuzi et al.l (j2003[ ).^ 
In short, we subtracted overscans and trimmed each im- 
age, and we combined and subtracted bias frames to re- 
move any residual bias structure. We also combined and 
applied dome flat fields in each filter to the science im- 
ages. We refined the image world coordinate systems 
using coordinates from the Guide Star Catalog 2 and 
projected the images onto the tangent plane, thereby 
resampling each pixel and correcting for the variable 
pixel scale. At this point in reduction we deviated from 
the NDWFS reduction techniques. In order to obtain 
the most accurate photometry and to acc ount correctly 
for ch ip-to-chip variations in color terms (jSlesnick et al.l 
|2002( ). we divided each resampled image into individual 
subimages for each chip; each subimagc was analyzed in- 
dependently. 

Object detection and photometric measurements were 
perfori ned using the D AOPHOT II suite of analysis pro- 
grams (jStetsonl Il987l ) and a detection threshold of 4(t. 
Approximately 50 bright, isolated stars were used in each 
subimage to define a point-spread function (PSF); pho- 
tometry for each star was determined via PSF-fitting. 
Curve-of-growth analysis of the PSF stars determined 
the aperture corrections used to convert PSF magnitudes 
to total instrumental magnitudes. 

The photometric catalogs from each band were then 
matched to create a final photometric catalog. Stars were 
required to have been detected on at least one frame in 
each of the three bandpasses in order to make the final 
catalog. Due to high background extinction, star-galaxy 
separation was deemed unnecessary. 

2.3. Calibration 

Calibrations for the MOSAIC data were determine d 
via imaging of standard star fields from iLandold (fl99l . 

IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
* Available at 



|http://www. noao.edu/noao/noaodeep/ReductionOpt/fr ames.html 
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V B-V U-B 

Fig . 1. — Difference between photometry of ISung &: Bessell 
1IT999I) and our photometry. Only slight systematics are observed 
in the B and \/-band photometry, while our U — B colors are sys- 
tematically bluer for hotter stars. 

The following transformation equations were used for 
calibration of the data: 



u = U + 2.5\ogt^^p + AQ + Ai{U-B) (1) 
+A2{X ~1.25) + A3{U-B)'^ 

b = B + 2.5logt,^p + Bo + Bi{B-V) (2) 
+B2{X - 1.25) 

v = V + 2.5\ogt^^p + Co + Ci{B-V) (3) 
+C2(X- 1.25), 

where u, b, v are the total instrumental magnitudes, 
U, B, V are the standard system magnitudes, tcxp is the 
exposure time, and X is the airmass. Examination of the 
standard star observation residuals revealed the necessity 
of the quadratic term in the ?7-band transformation. The 
number of observed standard stars with extreme colors 
was fairly small, rendering the value of the quadratic 
coefficient uncertain. However, the agreement of the 
WD [/-band photometry with theoretical values gives 
us confidence that this value is roughly correct; solv- 
ing the transformation equation without the quadratic 
term leads to observed WD colors in conflict with the 
evolutionary models. Values for the transformation co- 
efficients are given in Table [TJ 

As the number of observed standard stars was insuffi- 
cient to determine individual color terms for each CCD 
in the Mosaic camera, a single, mean value was deter- 
mined. For purposes of comparison, a second calibra- 
tion was performed a dopting published color t erms for 
each individual CCD (|Massev fc Slesnick|[l999l ): the re- 
sulting photometry for individual stars in B and V aver- 
aged 0.01 mag fainter with a dispersion of 0.01 mag. As 
iMassev fc Slesnidd (|1999l) did not use a quadratic term in 
the U -band transformati on, direct comparison of [/ -band 
terms is not possible, but iMassev fc Slesnic3 ()1999D note 
the chip-to-chip color-term variation introduces a scatter 
of ~ 0.04 mag over a 1 mag range in U — B. 

As a further sanity check, we compare our pho- 
tometry with previou s CCD photometry of M35 from 
ISung fc Besselll ()1999l ). shown in Figure [TJ We find no 
systematic offsets in the B or V photometry at > 0.02 
mag levels. However, a significant systematic trend is 
observable in U : for the bluest objects, we obtain U 
magnitudes systematically brighter by « 0.03 mag. It 
is unclear if this systematic is intrinsic to our data, per- 
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TABLE 1 

Photometric Transformation Equation Coefficients 



Filter Zero Point Color Term Airmass Term Quadratic Term 

U Ao = -22.981 ±0.009 Ai = -0.126 ± 0.021 A2 = 0.395 ± 0.059 A3 = 0.085 ± 0.019 

B Bo = -24.984 ± 0.006 Bi = -0.097 ± 0.005 B2 = 0.205 ± 0.031 

V Co = -25.157 ±0.004 Ci = 0.041 ± 0.004 C2 = 0.078 ± 0.021 




B-V U-B 

Fig. 2. — Color-magnitude diagrams for the field of M35. The 
thick, solid line is a 175-Myr, Z = 0.013 Padova isochrone. White 
dwarf cooling models are plotted for DA (solid) and DB (dotted) 
WD models with masses O.4Af0, O.8M0, and 1.2M0 (brighter to 
fainter) for cooling ages < 200Myr; DB models are only plotted for 
TcS < 30000K. The I.2M0 DB model is omitted, as it appears as 
a single point nearly coincident with the faint end of the 1.2Af0 
DA model. All models are shifted to the cluster distance and red- 
dening. White dwarf photometric selection criteria are shown as 
long-dashed lines; selected white dwarf candidates are indicated by 
squares. Error bars in white dwarf photometry are smaller than 
the point size. 

haps indicative that the color terms in Eq. [5] are a poor 
approximation to the tr ue color respo nse (see discussions 
oiU — B photom etry in 'Bessell' 'IPOS'), if the systematic 
is intrinsic to thelSung^ fc Bcsscll (1999) data, or both. 

2.4. Checking the Cluster Distance, Reddening and 
Metallicity 

The B — V,V and U — B,V color-magnitude diagrams 
(CMDs) of the stars in the field of M35 are shown in 
Figure [H Overplotted is the Padova isochrone for a 
Z = 0.01 3 ([Fe/H] = -0 .2 for the traditional solar abun- 
dances of I Anders fc Grev cssc 1989), 200-Myr-old stellar 
population shifted to the adopted cluster distance and 
reddening (see Hl.ip . The isochrone is found to be an 
excellent representation of the main sequence for V < 17 
(My < 6.5) in the B-V CMD. 

The same Z = 0.013, 200-Myr isochrone is too blue in 
the U-B CMD; agreement for V < 16 (My < 5.5) can 
be reached by increasing the reddening to E{B — V) = 
0.3. Such a change in reddening throws off the observed 
agreement in the B-V CMD. This suggests either a non- 
standard reddening law in the direction of M35, a signif- 
icant error in our [/-band zeropoint (in a direction which 
would worsen the systematic offset mentioned in the pre- 
vious section), or a problem with the isochrone U-hand 



calculations. Additional data and analysis are needed to 
explore this discrepancy further. 

A comparison of our data with Z = 0.008 and Z = 
0.019 Padova isochrones finds that the observed main se- 
quence agrees with the isochrone if the distance modulus 
of the cluster is shifted 0.15 mag closer or further, re- 
spectively. As with the Z — 0.013 isochrone, the U-B 
CMD requires a higher assumed reddening value to bring 
the isochrone and main sequence into agreement. 

The observed main sequence is also in excellent agree- 
ment with the e mpirical zero-age main sequence of 
iMermilliodI (|1981h . shifted to the adopted cluster dis- 
tance and reddening, and w i th th e fiducial M35 main 
sequence of Ivon Hippel et al.l (|2002D . 

We therefore conclude that our photometry is consis- 
tent with the literature-based distance, reddening, and 
metallicity we adopted in Hl.U We emphasize that, as 
our photometric data saturate well below the main se- 
quence turnoff, our main sequence data are not useful 
for age determination. 

2.5. Candidate WD Selection 

The CMDs spanning the entire magnitude range of 
our photometry are shown in Figure [2] along with cool- 
ing curves for hydrogen-atmosphere (DA) and helium- 
atmosphere (DB) WDs for a range of WD masses, all 
shifted to the adopted cluster distance modulus and red- 
dening. Several faint, blue objects are seen in the region 
of the CMD that should be populated by cluster WDs. 
The U-B, B-V color-color plot (Figure[3]) more clearly 
shows those objects with colors consistent with WD cool- 
ing models. 

Candidate WDs are selected using the criteria V > 18, 
B-V < 0.6, and U — B< —0.3. These criteria would 
include all WDs at the cluster distance and reddening 
except for the hottest low-mass (AfwD S O.4M0, Tcff > 
50,000K) WDs. These criteria are indicated in Figures 
[2] and [U candidate WDs selected via these criteria indi- 
cated as larger solid points in the CMDs and color-color 
plot. Coordinates, photometry, and cross- identifications 
for each candidate WD are given in Table [2l 

While most of the WD candidates are fairly isolated, 
some have close optical companions. In order to clarify 
identification for these WD candidates, finder charts for 
objects with close companions are provided in Figure |4l 

3. WD SPECTROSCOPY 

3.1. Observations and Data Reduction 

Spectra of several WD candidates were obtained be- 
tween 2002 December and 2005 November wit h the blue 
came r a of the LRIS spectrog raph on Keck I ()Oke et al.l 
119951 iMcCarthv erall 119981 ). We selected the 400 
grooves mm" 3400 A blaze grism for the highest avail- 
able throughput for the vital higher-order Balmer lines. 
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TABLE 2 

Candidate White Dwarfs in NGC 2168 



Object 




RA 






Dec 




V 


5V 


B - 


- V 


5{B 




U- 


-B 


S{U-B) 


LAWDS 


1 


6: 


:08: 


:38, 


.79 


24: 


:15: 


:06, 


.9 


20, 


.989 


0.019 


-0, 


.035 


0, 


.028 


-1, 


.095 


0, 


.028 


LAWDS 


2 


6: 


:08; 


:42, 


.30 


24: 


:10: 


:17, 


.7 


21 


.569 


0.032 


0. 


,061 


0, 


.044 


-1, 


.097 


0, 


.042 


LAWDS 


3 


6: 


:09; 


:04, 


.78 


24: 


:21: 


:39, 


.2 


20, 


.581 


0.020 


0. 


,053 


0, 


.028 


-0, 


.950 


0, 


.028 


LAWDS 


4 


6: 


:09: 


:05, 


.76 


24: 


:12: 


:11, 


.8 


21 


.216 


0.021 


0. 


148 


0, 


.030 


-0, 


.858 


0. 


.030 


LAWDS 


5 


6: 


:09: 


:11, 


.54 


24: 


:27: 


:20, 


.9 


20, 


.065 


0.017 


-0, 


.128 


0, 


.024 


-1, 


.173 


0. 


.025 


LAWDS 


6 


6: 


:09: 


:23, 


.48 


24: 


:27: 


:22, 


.0 


19, 


.863 


0.016 


-0, 


.128 


0, 


.023 


-1, 


.167 


0, 


.024 


LAWDS 


10 


6: 


:09: 


:43, 


.63 


24: 


:19: 


:15, 


.8 


20, 


.238 


0.020 


0. 


122 


0, 


.030 


"0, 


.506 


0, 


.032 


LAWDS 


11 


6: 


:09: 


:42, 


.79 


24: 


:11: 


:05, 


.4 


21 


.198 


0.025 


0. 


110 


0, 


.037 


-0, 


.680 


0. 


.038 


LAWDS 


12 


6: 


:09: 


:31, 


.19 


24: 


:19: 


:06, 


.2 


21 


.303 


0.026 


0. 


,045 


0, 


.038 


-1, 


.009 


0, 


.038 


LAWDS 


13 


6: 


:09: 


:29, 


.71 


24: 


:15: 


:58, 


.6 


21 


.542 


0.024 


0. 


,341 


0, 


.037 


-0, 


.551 


0, 


.038 


LAWDS 


14 


6: 


:09: 


:15, 


.10 


24: 


:33: 


:15, 


.4 


21 


.701 


0.027 


0. 


,059 


0, 


.038 


-0, 


.933 


0. 


.036 


LAWDS 


15 


6: 


:09: 


:11, 


.63 


24: 


:02: 


:38, 


.5 


20, 


.785 


0.022 


-0, 


.039 


0, 


.032 


-1, 


.019 


0, 


.033 


LAWDS 


18 


6: 


:08: 


:35, 


.44 


24: 


:34: 


:19, 


.5 


22, 


.566 


0.043 


0. 


172 


0, 


.071 


-0, 


.671 


0, 


.069 


LAWDS 


22 


6: 


:08: 


:24, 


.65 


24: 


:33: 


:47, 


.6 


19, 


.657 


0.016 


0. 


,008 


0, 


.023 


-0, 


.944 


0, 


.025 


LAWDS 


27 


6: 


:09: 


:06, 


.26 


24: 


:19: 


:25, 


.3 


21 


.398 


0.026 


0. 


,090 


0, 


.039 


-1, 


.009 


0, 


.039 


LAWDS 


28 


6: 


:08; 


:13, 


.50 


24: 


:20: 


:32, 


.5 


21 


.631 


0.029 


0. 


,018 


0, 


.042 


-1, 


.370 


0, 


.040 


LAWDS 


29 


6: 


:08; 


:02, 


.20 


24: 


:25: 


:24, 


.2 


20, 


.719 


0.019 


-0, 


.075 


0, 


.027 


-1, 


.087 


0, 


.027 


LAWDS 


30 


6: 


:07: 


:56, 


.63 


24: 


:13: 


:27, 


.2 


21 


.175 


0.020 


-0, 


.005 


0, 


.029 


-0, 


.927 


0. 


.030 


LAWDS 


31 


6: 


:10: 


:08, 


.48 


24: 


:22: 


:32, 


.3 


20, 


.231 


0.020 


0. 


,048 


0, 


.031 


-0, 


.991 


0. 


.033 


LAWDS 


32 


6: 


:09: 


:37, 


.35 


24: 


:31: 


:52, 


.5 


22, 


.000 


0.030 


0. 


,418 


0, 


.047 


-0 


.432 


0, 


.062 


LAWDS 


33 


6: 


:09: 


:33, 


.02 


24: 


:15: 


:23, 


.5 


22, 


.259 


0.033 


0. 


,284 


0, 


.051 


-0, 


.488 


0. 


.060 


LAWDS 


34 


6: 


:09: 


:25, 


.26 


24: 


:14: 


:05, 


.2 


21 


.026 


0.022 


0. 


,550 


0, 


.033 


-0, 


.495 


0, 


.035 


LAWDS 


35 


6: 


:09: 


:12, 


.92 


24: 


:22: 


:17, 


.4 


21 


.809 


0.027 


0. 


,342 


0, 


.040 


-0 


.402 


0, 


.042 


LAWDS 


36 


6: 


:09: 


:04, 


.59 


24: 


:06: 


:45, 


.7 


21 


.276 


0.021 


0. 


,465 


0, 


.033 


-0, 


.483 


0. 


.036 


LAWDS 


37 


6: 


:08; 


:59, 


.02 


24: 


:08: 


:40, 


.8 


21 


.049 


0.020 


0. 


,406 


0, 


.030 


-0 


.554 


0, 


.033 


LAWDS 


38 


6: 


:09; 


:08, 


.26 


24: 


:36: 


:24, 


.2 


22, 


.158 


0.034 


0. 


,277 


0, 


.051 


-0 


.444 


0, 


.056 


LAWDS 


39 


6: 


:09: 


:01, 


.46 


24: 


:26: 


:50, 


.3 


22, 


.648 


0.052 


0. 


,275 


0, 


.081 


-0 


.790 


0, 


.086 


LAWDS 


40 


6: 


:08: 


:54, 


.53 


24: 


:35: 


:58, 


.5 


22, 


.244 


0.043 


0. 


,222 


0, 


.061 


-0, 


.612 


0. 


.070 


LAWDS 


41 


6: 


:08; 


:34, 


.97 


24: 


:32: 


:47, 


.3 


18, 


.886 


0.015 


0. 


,471 


0, 


.022 


-0 


.473 


0, 


.024 


LAWDS 


42 


6: 


:08; 


:24, 


.11 


24: 


:22: 


:34, 


.7 


22, 


.338 


0.042 


0. 


,277 


0, 


.069 


-0, 


.484 


0, 


.076 


LAWDS 


43 


6: 


:08; 


:20, 


.96 


24: 


:08: 


:51, 


.4 


22, 


.440 


0.037 


0. 


,523 


0, 


.065 


-0, 


.456 


0, 


.086 


LAWDS 


44 


6: 


:08: 


:03, 


.84 


24: 


:27: 


:37, 


.2 


19, 


.718 


0.017 


0. 


167 


0, 


.025 


-0, 


.433 


0, 


.027 


LAWDS 


45 


6: 


:08; 


:03, 


.77 


24: 


:27: 


:38, 


.3 


18, 


.403 


0.015 


-0, 


.122 


0, 


.021 


-0, 


.929 


0, 


.022 


LAWDS 


46 


6: 


:08: 


:00, 


.63 


24: 


:07: 


:40, 


.2 


20, 


.998 


0.019 


0. 


176 


0, 


.029 


-0, 


.575 


0. 


.030 


LAWDS 


47 


6: 


:08: 


:00, 


.37 


24: 


:18: 


:02, 


.1 


21, 


.540 


0.023 


0. 


,422 


0, 


.036 


-0, 


.521 


0. 


.038 


LAWDS 


48 


6: 


:07; 


:58, 


.70 


24: 


:28: 


:40, 


.1 


21, 


.590 


0.030 


0. 


,336 


0, 


.047 


-0, 


.567 


0, 


.048 


LAWDS 


49 


6: 


:07: 


:52, 


.95 


24: 


:25: 


:23, 


.5 


20, 


.191 


0.017 


0. 


133 


0, 


.025 


-0, 


.970 


0. 


.026 


LAWDS 


50 


6: 


:07: 


:47, 


.77 


24: 


:28: 


:52, 


.0 


21, 


.465 


0.030 


0. 


,300 


0, 


.044 


-0 


.477 


0. 


.046 


LAWDS 


51 


6: 


:07; 


:47, 


.39 


24: 


:34: 


:15, 


.3 


21, 


.877 


0.031 


0. 


182 


0, 


.044 


-0 


.680 


0, 


.041 


LAWDS 


52 


6: 


:07: 


:47, 


.37 


24: 


:35: 


:23, 


.5 


20, 


.653 


0.018 


0. 


140 


0, 


.027 


-0 


.764 


0, 


.028 


LAWDS 


53 


6: 


:07; 


:33, 


.79 


24: 


:07: 


:55, 


.5 


21, 


.341 


0.028 


0. 


,560 


0, 


.046 


-0, 


.442 


0, 


.050 



Cross-Identification 



NGC 2168-1 in Romanishin fc Angel (1980) 

NGC 2168-2 in Romanishin & Angel (1980) 
NGC 2168-3 in Romanishin & Angel (1980) 
NGC 2168-4 in Romanishin fc Angel (1980) 



Discovered by ^oi^jJi££e^_et_^l^ (2002 



V57 in Mocheiska et al. (2004 



Note. — Units of right ascension are hours, minutes and seconds, am 
are for Equinox J2000.0. The official (i.e., lAU-approved) format for e 

The D560 dichroic was used to obtain simultaneous ob- 
servations of the Ha Hne (not presented here) . Our aper- 
ture was a 1" longsht at paraUactic angle. The resulting 
spectral resolution (FWHM) is « 6A. 

We reduced the spectra using the onedspec package in 
IRAF. Overscan regions were used to subtract the ampli- 
fier bias. Flat fielding was complicated by the discovery 
of a sharp inflection point in the response at « 4200A 
and two low-level (w 3%) emission lines in the illumi- 
nated dome flat field spectra between 3950A and 3975A; 
both of these features introduced ringing into the flat 
field. The ringing due to the inflection point was elimi- 
nated by creating a piecewise-smooth response function 
for image sections below and above the inflection point. 
As internal flat fields lacked the emission features, only 
internal flats were used to create the final fiat fields. 

Cosmic rays were removed from the two-dimensional 
spectrum using th e "L.A. Cosmic" La placian cosmic-ray 
rejection routine (jvan DokkumI [200ll ) . We then aver- 
aged multiple exposures of individual objects weighted by 
their individual signal-to-noise ratios and extracted the 
one-dimensional spectrum. We applied a wavelength so- 
lution derived from Hg, Cd, and Zn arclamp spectra. We 



units of declination are degrees, arcminutes, and arcseconds. Coordinates 
,ch object name is NGC 2168:LAWDS NN. 

determined and applied a relative flux calibration from 
long-slit spectra of multiple spectrophotometric standard 
stars. We made no attempt at obtaining absolute spec- 
trophotometry for any object. 

The observed candidates, their dates of observation, 
and the spectroscopic identification of each object are 
given in Table O all spectra are shown in Figure [5l 
NGC 2168:LAW DS 28 is the hot DQ WD discussed in 
IWilliams et all (,200& ): NG C 2168 :LAWDS 4 is spectral 
type DB and is discussed in ij4.6.2l All twelve other WDs 
are of spectral type DA. 

Figure [5] shows the WD regions of the CMDs and color- 
color plot of M35 indicating each object's spectroscopic 
identification. In the B — V vs. V CMD, the cluster 
WDs form a very tight sequence about the IM© DA 
cooling model, with the exception of the likely magnetic 
DA LAWDS 22 (see next section). The redder color of 
LAWDS 22 could be due to an unresolved low- mass com- 
panion; the red side of LRIS was not functioning on the 
night this WD was observed, so we are unable to examine 
the spectra for evidence of such a companion. The hot 
DQ LAWDS 28 is significantly bluer inU-B than aU the 
other WDs. This exceptional color is due to high carbon 
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3 




B-V 

Fig. 3. — Color-color plot for the field of M35. Small points 
are for all detected objects; larger squares are the selected WD 
candidates. Cooling models are shown for t^ooi < 200Myr for DA 
WDs (solid) of masses OAMq, O.SMq, and 1.2Mq (bottom to top) 
and for DB WDs of masses OAMq and O.SM©; the 1.2Mq model is 
omitted. Models are shifted to the assumed cluster reddening. The 
dashed line indicates WD selection criteria. The arrow indicates 
the reddening vector of E{B — V) = 0.22, assuming Ry = 3.1. 



TABLE 3 

Spectroscopically-Observed Objects 



Object 


UT Obs Date 


Total exp. (s) 


S/N'' 


ID 


LAWDS 1 


2002 Dec 8 


2700 


64 


DA 


LAWDS 2 


2002 Dec 8 


2700 


109 


DA 




2004 Oct 11 


3600 






LAWDS 4 


2002 Dec 8 


2700 


32 


DB 




2004 Oct 11 


900 






LAWDS 5 


2002 Dec 8 


2700 


215 


DA 




2002 Dec 9 


3600 






LAWDS 6 


2002 Dec 8 


2700 


247 


DA 




2002 Dec 9 


2700 






LAWDS 11 


2004 Feb 12 


3600 


57 


DA 




2005 Nov 26 


1800 








2005 Nov 27 


4800 






LAWDS 12 


2005 Nov 26 


4800 


91 


DA 


LAWDS 14 


2005 Nov 26 


6000 


73 


DA 


LAWDS 15 


2004 Feb 12 


3300 


86 


DA 


LAWDS 22 


2004 Feb 12 


3300 


198 


DA(H?) 


LAWDS 27 


2002 Dec 9 


3600 


96 


DA 




2004 Feb 12 


3600 








2004 Oct 11 


3300 






LAWDS 28 


2004 Oct 12 


1200 


73 


DQ 




2005 Nov 26 


4800 






LAWDS 29 


2004 Oct 12 


3600 


75 


DA 


LAWDS 30 


2005 Nov 27 


6000 


41 


DA 



^ Average of signal-to-noise per resolution element at pseudocon- 
tinuum surrounding KS (for DA WDs) or at 4200A (non-DAs) 

opacity at wavelengths < 1500A, which re-distributes the 
ultraviolet flux into the near-UV (jPufour et al.ll2008[ ). 

3.2. WD parameter determination 

Toff and log 5 were determined for each DA WD via 
simu ltaneous fitting of th e H/3 to H9 Balmer line pro- 
files (jBergeron et al.iri992[ ). Model atmospheres used for 
this fitting are derived from Koester's model grids; de- 
tails of the in p ut ph ysics and methods can be found in 
iKoester et al.l (|2001[ ) and references therein. The algo- 



rithm fo r our fitting routin e and e rror estimation are de- 
tailed in IWilliams fc Boltd ()2007n . with the change that 
the model grid was expanded to include higher surface 
gravities (7.0 < \ogg < 10.0). Signal-to-noise was cal- 
culated by determining the root mean square scatter per 
GA resolution element about the fit pseudo-continuum on 
both sides of the H(5 line. 

Initial spectral fits proved unsatisfactory for the 
hottest white dwarfs (LAWDS 5, LAWDS 6, and LAWDS 
22), with the observed spectrum deviating strongly from 
the best-fit models. The models severely underestimated 
the depth of the He line and overestimated those of lower- 
order Balmer lines. There are several potential sources 
for the poor He fits in these three WDs: 

• Magnetic fields — In the spectrum of LAWDS 
22, the cores of the Balmer lines are signifi- 
cantly fiattened. This is likely to be due to 
a magnetic field with a strength sufficient to 
broaden the line without producing obvious Zee- 
man splitting at our spectral re solution ('^ 500 
kG; IWickramasinghe fc Ferrariol i2000( ): a rough 
comparison to a magnetic atmospheric models re- 
veals good agreement for a magnetic field strength 
of ~ 1 MG (A. Kanaan, private communication). 
The Balmer line models we fit do not include mag- 
netic fields, and so the parameters derived for this 
star are invalid. Weaker fields could also be present 
in LAWDS 5 and LAWDS 6, accounting for their 
poor fits. However, there is no compelling evidence 
for magnetism in these two stars. 

• Known high Toff systematics — DA WDs with 
TcS ^ 40000 K show systematic differences 
in Tcff derived from different methods, though 
the proposed causes of these differences, includ- 
ing non-LTE effects, small quantities of helium, 
and meta l opacities, remain a matter o f de- 
bate (e.g., iBer gcron et al. 1994; Naoiwotzki" '19971 : 
iBarstow et aLli2003;Taioie fc Bergeron 2007.) . The 
sense and magnitude of our deviations are s imi- 
lar to those calculated bv IBarstow et al.] ()1998[ ) for 
a metal-rich, non-LTE atmosphere compared to a 
pure-H, LTE model. 

• Interstellar absorption — In our highest signal-to- 
noise spectra, we observe an additional absorption 
feature on the red wing of the H7 line that we iden- 
tify as the 4430 A diffuse interstellar band (DIB; 
e.g., Hcrbig 1966). Weak Ca II K absorption is also 
visible in the higher signal-to-noise WD spectra; 
due to the high Toff of these WDs, this absorption 
must be interstellar. The He lines will therefore be 
contaminated by Ca II H absorption. In the hottest 
WDs, the depth of Ca II hues are - 50% of the He 
line. We have not attempted to model and sub- 
tract the Ca II H line, so this line will have some 
effect on the model atmosphere fits. However, as 
the Ca II H line is significantly narrower than the 
He line, and the deviations span the entire He line, 
contamination from Ca II H is therefore not the 
primary cause of the poor fits. 

• Shortcomings in the model atmospheres — The 
high-order Balmer lines are strongly affected by 
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Fig. 4. — Finder charts for white dwarf candidates with close optical companions, 
is down. Note that in panel (f), each component in an optical binary is identified. 



Each image is 1' on a side. North is to the left; east 
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LAWDS 28 



LAWDS 12 . 




^ LAWDS 1 



LAWDS 27 




_L 



3500 



4000 



4500 



5000 



A(A) 



Fig. 5.— Spectra of observed WDs in the field of NGC 2168. 
Spectra have been normalized at 5200A and arbitrary vertical off- 
sets applied. The top two spectra are hot DQ (LAWDS 28) and 
the DB (LAWDS 4); the remaining are DAs, ordered from top to 
bottom in decreasing Tefj. 
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23 
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Fig. 6. — The white dwarf regions of the M35 B-Y (left) and 
U — B (right) CMDs with spectroscopic identifications. Large filled 
dots are DA white dwarfs; the open circle is the hot DQ white 
dwarf, and the DB white dwarf is plotted as a diamond. Small 
points are white dwarf candidates not observed spectroscopically. 
Curves are cooling models for DA (solid) and DB (dashed) white 
dwarfs with cooling ages < 200Myr for masses of 0.4AfQ, O.GMq, 
O.SMq, I.OMq, and 1.2A/0 WDs (top to bottom). Curves are 
shifted to the white-dwarf derived distance and reddening values 
(see i|4.5l l. In the B-V CMD, the cluster WDs form a tight cool- 
ing sequence near the IMq cooling track, though these are more 
scattered in the U-B CMD. 



the dissolution of the higher atomic energy lev- 
els due to interactions with perturbing particles. 
This is described in the atmospheric models by 
the Hummer-Mihalas-Deppen occupation probabil- 
ity, which has free parameters that are difficult to 
quantify. The perturbation effect may be a bit too 
strong, leading to weaker higher-order lines in the 
models. To test this possibility, we have used our 
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models and fitting routine to determine parameters 
of hot, high- gravity WDs in the Palomar Green 
WD sample (iLiebert et al.ll20d5al ). whose spectra 
were graciously provided by J. Liebert. These fits 
do not show the same problems with the fitting 
of He observed here; therefore, we conclude that 
shortcomings in the models are likely not responsi- 
ble for our poor spectral fits. 

• Data reduction systematics — As mentioned in 
the dome spectroscopic flat field lamps were 
observed to have weak emission features in this 
wavelength region, so only internal flat field lamps 
were used. If the internal lamps have similar spec- 
tral features, the He profiles could be affected. 
However, no such features were apparent in the in- 
ternal flats, and the cooler DAs in this study do 
not show systematics in the He line, so we consider 
this explanation unlikely. 

In short, we find no compelling explanation why the He 
lines in LAWDS 5 and LAWDS 6 are not well-represented 
by the atmospheric models; the failure in LAWDS 22 
seems likely to be due to magnetic fields. 

As a check against model atmosphere or data reduction 
systematics affecting all of our WDs, we re-fit each WD 
excluding the He line from the fits. In all cases, the 
Teff values are identical within the internal fitting errors 
(described below). This is expected, as the majority of 
leverage in T^g determination is from the lower Balmer 
lines. The logg values are systematically higher by an 
average of 0.1 dex when the He line is excluded. This 
change is because the higher-order Balmer lines are most 
sensitive to surface gravity, and the absence of the H8 and 
H9 lines due to the dissolution of these energy levels in 
these higher gravity atmospheres only sets a lower limit 
on the surface gravity. 

The best-fitting model atmospheres from each of the 
two fits were plotted over each WD's observed spectrum, 
and the best qualitative fit was adopted for the remain- 
ing analysis; adopted values are given in Table ID For 
most WDs, this adopted fit was the fit obtained using all 
Balmer lines. For LAWDS 5, LAWDS 12, and LAWDS 
22, the adopted fit excluded He. For LAWDS 6, both 
fits were qualitatively similar, so an intermediate value 
(also qualitatively similar) was adopted, and the internal 
error bars were increased to include both fit values. The 
adopted-fit model Balmer lines are shown plotted over 
the observed spectra in Figure [71 

Masses (Mwd) and coohng ages (tcooi) for each white 
dwarf were computed from evolutionary models pro- 
vided by P. Bergeron & G. Fontaine. These models 
assume "thick" hydrogen layers {Mh — 10^^Af„ ). For 
Tcff > 30000K, the pure carbon cooling models of lWoodI 
( 11995|| are used; for T e g < 3 0000K, the mixed C/0 mod- 
els of iFontaine et al.] ()2001[ ) are used. The evolutionary 
models al so provide synthet ic photometry as described in 
iHolberg fc BergeronI ()2006f ). As pure carbon models are 
generally disfavored at these masses, we also calculated 
Mwd and TcooI f or T^m > 300 00K WDs using the C/0 
mixed models of IWoodI ( 1995 ?): the resulting differences 
in Mwd and TcooI were within the stated error bars. The 
derived Mwd and TcooI for each WD are presented in Ta- 
ble H Due to the poor fits of LAWDS 5 and LAWDS 6, 
their derived masses are of questionable quality. 




-100 100 -100 100 -100 100 

AX 

Fig. 7. — Balmer line profiles and best-fit models for 12 DA white 
dwarfs in the field of M35, from H/3 (bottom) to H9 (top). Plotted 
models are the adopted best-fit models from Table [3] Significant 
deviations in He from the best-fit models are noticeable in LAWDS 
6, LAWDS 5, and LAWDS 22. The cores of the fines for LAWDS 22 
are noticeable flattened, likely indicative of nearly-resolved Zeeman 
splitting. The 4430A diffuse interstellar band can be seen in the 
red (right) wing of H7 in many of the higher signal-to-noise WDs, 
and interstellar Ca II K can be often seen in the blue wing of He. 
WDs are arranged in order of decreasing T^g. 



3.2.1. Error Determinations 

We consider two primary sources of error in the Tog 
and logg values determined for each WD: 

Internal fitting errors — We refer to the random er- 
rors in our determination of Tcs and log g resulting from 
the signal-to-noise of our observations as "internal fit- 
ting errors," or the distribution of spectral parameters 
we would expect given an ensemble of observations with 
the same signal-to-noise of the same WD using the same 
instrumental setup. We determine the magnitude of 
these errors via th e Mon te Carlo technique described in 
iWilUams &: Boiti ((2007h . For LAWDS 5, LAWDS 12 
and LAWDS 22, where the adopted atmospheric param- 
eters used the fits excluding the He line, that line was 
also excluded in the Monte Carlo error determination, 
thereby providing a realistic estimate of the additional 
random error introduced by the exclusion. 

External fitting errors — It has been noticed that 
the scatter in measured Toff and logg derived by dif- 
ferent groups for the same WDs is larger than the in- 
ternal errors calculated by Balm er fitting routin es (e.g., 
iNapiwotzki et al.lll999t lWimams~fc Boltel 120071 ). Possi- 
ble sources for these differences include the different in- 
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TABLE 4 

DA White Dwarf Spectral Fits and Derived Parameters 



Object 




log 9 


Internal Errors 


Total Errors 




<5MwD 


T-cool 


cool 


Minit 




(K) 






<51ogg 




<51ogg 




(^%) 


log(yr) 


log(yr) 




IjAW Do 1 


33400 


8.36 


320 


0.11 


1150 


0.16 


0.873 


0.091 


7.228 


0.260 


4.39 




Q/1 1 nn 
o4iUU 


O.DZ 


OIU 


n OA 


1140 


U. lO 


i.Uio 


U.UD i 


/ .00 i 




4. / y 


LAWDS 5^^ 


SSTSO"* 


8.39=* 


850 


0.05 


1390 


0.13 


0.916 


0.075 


6.103 


0.088 


4.21 


LAWDS 6=^ 


56000== 


8.32° 


1500 


0.14 


1860 


0.18 


0.877 


0.096 


6.077 


0.082 


4.21 


LAWDS 11 


20800 


8.28 


520 


0.08 


1260 


0.16 


0.802 


0.096 


8.047 


0.173 


6.63 


LAWDS 12 


34500=* 


8.44=* 


280 


0.10 


1140 


0.16 


0.922 


0.092 


7.314 


0.281 


4.43 


LAWDS 14 


29200 


8.62 


220 


0.07 


1120 


0.14 


1.010 


0.072 


7.865 


0.167 


5.32 


LAWDS 15= 


29500 


8.40 


160 


0.00 


1110 


0.15 


0.888 


0.088 


7.551 


0.236 


4.63 


LAWDS 22" 


50000=* 


8.08=* 


1420 


0.07 


1800 


0.12 




f 


f 


f 


f 


LAWDS 27 


30200 


8.64 


100 


0.04 


1100 


0.13 


1.022 


0.072 


7.840 


0.170 


5.22 


LAWDS 29 


32400 


8.38 


220 


0.07 


1120 


0.14 


0.882 


0.078 


7.331 


0.235 


4.44 


LAWDS 30 


30400 


8.62 


390 


0.18 


1170 


0.22 


1.011 


0.122 


7.808 


0.289 


5.12 



(5Mi 



'imt.Obs 


oJWinit.Sys 






+0.23 


+0.35 


-0.09 


-0.27 


+0.47 


+0.46 


-0.26 


-0.36 


+0.00 


+0.29 


—0.00 


—0.23 


+0.00 


+0.29 


-0.00 


-0.23 


+0.00 


+ 1.89 


-1.38 


-0.96 


+0.30 


+0.36 


-0.13 


-0.27 


+0.95 


+0.73 


-0.44 


-0.48 


+0.38 


+0.41 


-0.22 


-0.32 


f 


f 


+0.71 


+0.68 


-0.42 


-0.45 


+0.23 


+0.36 


-0.12 


-0.28 


+ 1.49 


+0.63 


-0.55 


-0.43 



^ Error due to total fitting crror*^ Additional systematic error due to cluster age uncertainty^^ Re-reduction of data presented in IWilliams et all 
|20p4)=* Fit excludes H€== Adopted values intermediate to fits obtained by including and excluding He*^ Values not calculated due to poor Balmer 
line fits 



struments with which these stars were observed, the dif- 
ferent spectral fitting routines, and the differen t model 
atmospheres used. In IWilliams fc Bolti ()2007D . we de- 
termined that these external errors appear to be random 
with a scatter of « HOOK in Toff and a scatter of « 0.12 
dex in logg based on a comparison of our observations 
of field DA WD spectra to fits i n the literature. In our 
past work (with the exception of lRubin et al.l[2008[ ). we 
have not propagated the external fitting errors through 
our analysis. In order to facilitate proper comparison 
of our open cluster WDs with those of other groups, we 
propagate these errors through the analysis in this paper. 

We obtain the total fitting error in Tea and log g for 
each WD by adding the internal, and external fitting er- 
rors in quadrature. This total error is then propagated 
through the analysis to obtain the stated errors in all sub- 
sequently derived parameters (Mwd, t'cooI, and Minit)- 
The internal and total fitting errors are both given in 
Table m 

4. THE WD POPULATION OF M35 
4.1. Cluster membership 

In the absence of kinematic information, such as radial 
velocity or proper motion measurements, the best means 
of determining cluster membership of WDs is to apply 
age and distance criteria. Any current cluster member 
WDs must have TcooI shorter than the cluster age and, 
unless escaped from the cluster, must lie at the same 
distance. Both TcooI and My for each WD are deter- 
mined from the evolutionary models. For each WD, TcooI 
is shown (along with other derived mass quantities) in 
Table m and (m — M)y is given along with other derived 
photometric quantities in Table El We identify candi- 
date cluster member WDs as those DAs with TcooI < 200 
Myr and that have distance moduli consistent within 2cr 
of the cluster distance modulus. All of the WDs in this 
sam ple meet these c r iteria . 

In IWilliams eTall ((200l . we suggested that LAWDS 
15 could also be an escaped cluster member. After the re- 
analysis presented in this paper, the WD distance mod- 
ulus is consistent with the cluster distance modulus, so 
the speculation on its potential escape from the cluster 



was premature. 

There are some factors that lead us to question the 
membership of DA LAWDS 11. Its cooling age is ~ 
40 Myr older than the second-oldest WD, LAWDS 14, 
though this difference is within the stated errors. In ad- 
dition, LAWDS 11 is the second-least massive WD in 
the cluster, despite having the highest progenitor mass. 
This goes against the preconceived notion that higher- 
mass stars produce higher-mass WDs. There are theo- 
retical means to explain this low mass, such as binary 
evolution or enhanced mass loss during the post-main 
sequence evolution of the star. 

It is also possible that LAWDS 11 is a field WD un- 
related to the star cluster that lies within our selection 
v olume. U sing a calculation similar to that described in 
R ubin et a l. (2008), we estimate « 3.2 field WDs meet- 
ing our cluster-member criteria should be in our imaged 
area. As ~ 22% of field WDs in a, volume-limit ed sam- 
ple are more massive than CSM^ llLiebert et al.| [2C)05a). 
it would not be improbable to find a massive field WD 
in our sample. Measurements of this star's proper mo- 
tion are needed to determine if it is indeed a member 
of the cluster; as LAWDS 11 meets all of our current 
membership criteria, we include it in our analysis below. 

4.1.1. Sample Completeness 

The fourteen WDs observed in this cluster may not be 
the complete WD population of M35. Our imaging is 
only of the central w 18 arcmin radius o f the cluster; the 
tidal radius of M35 is > 33 arcmin (jLeonard fc MerrittI 
I1989D . Therefore, it is fully possible that cluster WDs 
may be found outside our imaged area. However, al- 
most all main-sequence proper motion members more 
massive than 1.2M0 are l ocated within a 20 arcmin ra- 
dius of the cluster center (iMcNaniara fc Sekiguchilll986l : 
iBarrado v Navascues et al.l l2001bl ). and the progenitor 
stars of M35 WDs were significantly more massive than 
this, suggesting that few cluster member WDs may be 
outside of our images. 

We may also be missing cluster WDs within our im- 
aged area. Because precise WD parameter measure- 
ments rely on the higher-order Balmer lines and because 
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TABLE 5 

Model-dependent photometric properties of NGC 2168 white dwarfs 



Model Photometry Derived Distance and Reddening 



WD 


Mv 






"i.B-V)o 


{U-V)o 


"(U-V)o 


{m~M)i/ 


^{m-M)v 


E{B-V) 


''e(B-V) 


E{U-V) 


''E{U-V) 


1 


10.300 


0.283 


-0.221 


0.011 


-1.420 


0.016 


10.689 


0.284 


0.186 


0.030 


0.290 


0.043 


2 


10.736 


0.254 


-0.222 


0.010 


-1.433 


0.014 


10.833 


0.256 


0.283 


0.045 


0.397 


0.062 


5 


9.759 


0.247 


-0.291 


0.003 


-1.563 


0.004 


10.306 


0.248 


0.163 


0.024 


0.262 


0.035 


6 


9.592 


0.319 


-0.295 


0.003 


-1.570 


0.004 


10.271 


0.319 


0.167 


0.023 


0.275 


0.034 


11 


11.090 


0.270 


-0.051 


0.025 


-1.006 


0.045 


10.108 


0.271 


0.161 


0.045 


0.436 


0.069 


12 


10.382 


0.300 


-0.228 


0.009 


-1.438 


0.013 


10.921 


0.301 


0.273 


0.039 


0.474 


0.055 


14 


11.041 


0.264 


-0.169 


0.016 


-1.324 


0.025 


10.660 


0.265 


0.228 


0.041 


0.450 


0.058 


15 


10.622 


0.270 


-0.179 


0.015 


-1.330 


0.025 


10.163 


0.271 


0.140 


0.035 


0.272 


0.052 


27 


11.008 


0.265 


-0.182 


0.014 


-1.353 


0.021 


10.390 


0.266 


0.272 


0.041 


0.434 


0.059 


29 


10.393 


0.246 


-0.212 


0.011 


-1.401 


0.017 


10.326 


0.247 


0.137 


0.029 


0.239 


0.042 


30 


10.957 


0.428 


-0.185 


0.016 


-1.358 


0.024 


10.218 


0.428 


0.180 


0.033 


0.426 


0.048 



an atmospheric dispersion corrector was not available 
when these observations were made, we decided to make 
longslit observations of individual WDs at parallactic an- 
gle rather than multi-slit observations at less favorable 
angles. Therefore, most of the WD candidates lack spec- 
troscopic observations. 

Many of the unobserved WD candidates lie along and 
above the 0.6Mq WD cooling curve; all of the observed 
DAs are more massive than this and, with the excep- 
tion of LAWDS 22 (see H3.1\i . lie along more massive 
WD cooling curves. We cannot rule out that the unob- 
served candidates may be cluster WDs in double degener- 
ate or WD-I-M binaries, or low-mass cluster WDs formed 
via some binary process. While interesting in their own 
right, any of these scenarios would warrant their exclu- 
sion from the initial-final mass relation and A/up analyses 
below, which are assumed to be valid only for single-star 
evolution. 

The other unobserved WD candidates, if cluster mem- 
bers, would have cooling ages > 200 Myr, the oldest 
likely age for the cluster. Therefore, they are almost cer- 
tainly not white dwarfs related to the star cluster. 

Most of the WD candidates lacking spectroscopic ob- 
servations are likely field WDs. As mentioned in the 
previous section, we estimate that, on average, 3.2 field 
WDs meeting our cluster membership selection criteria 
should be found within our field of view. Our photomet- 
ric selection criteria are much more generous, consistent 
with WDs as cool as ~ 7000K at the cluster distance 
and reddening; such WDs would be significantly older 
than the star cluster, and therefore not cluster members. 
Solely using our photometric selection criteria and inte- 
grating along the line of sight to the cluster, we estimate 
~ 35 field WDs to be in our field of view. Given our 
41 total WD candidates, 10 to 14 of which are cluster 
members, we expect that most or all of these unobserved 
WD candidates are field WDs. 

In the B-V vs. V CMD (Fig. El), a tight, nearly con- 
tinuous sequence of member WDs is observed. There 
are no remaining WD candidates along the IMq cooling 
track younger than 200 Myr, suggesting that no older, 
massive WDs are present in the cluster. However, we 
cannot rule out that cooler WDs were formed and lost 
from the cluster due to dynamical processes. Photomet- 
ric incompleteness is also a concern; artificial star tests 
(significantly complicated by proper treatment of the mo- 
saicked CCDs) are being performed and will be presented 
in a future paper. 



In short, spectroscopic observations of the unobserved 
WD candidates are needed to ascertain their true nature, 
but they are quite unlikely to impact the conclusions of 
this paper. 

4.2. The initial-final mass relation 

The progenitor star masses for cluster WDs can be de- 
termined by subtracting TcooI from the cluster age. The 
result represents the lifetime of the progenitor star from 
the zero-age main sequence through the planetary nebula 
stage. We then use the jFe/H| = —0.2 evolutionary mod- 
els of iMarigo fc GirardU (|2007D to determine the mass of 
star with this lifetime. These model lifetimes only ex- 
tend from the zero-age main sequence through the start 
of the thermally-pulsing AGB phase. The amount of 
time from the first thermal pulse through the thermally- 
pulsing AGB and planetary nebula phases to the start of 
the WD cooling track is negligible (~ 10^ yr) compared 
with the total stellar lifetimes 10^ yr) and the errors 
on the cluster age 10^ yr). 

Table m gives the progenitor mass Minit for each cluster 
WD. This table excludes parameters for LAWDS 22 due 
to its likely magnetic nature and resulting poor spectral 
fits. The table lists the observational errors in Mjnit re- 
sulting from the propagation of the total fitting errors 
through the calculations. Also given are the systematic 
errors, the changes in Mjnit due to the 25 Myr uncertainty 
in the cluster age. These systematic errors would be in 
the same sense for each cluster WD and so should not be 
added in to the random error (see ij4.3p . The given errors 
in Minit are strongly asymmetric and represent the 68% 
confidence level. The initial-final mass relation is plot- 
ted in Figure [5] Error bars are the 68% confidence levels 
derived from the total fitting errors. Despite the uncer- 
tainty in the spectral fits for LAWDS 5 and LAWDS 6, 
their Minit values and uncertainties are fairly robust; as 
their TcooI are very short, even an error in Toff of 10000 
K translates to an error in Minit of only 0.02Mq. 

Also included in Figure [8] are cluster WDs from the 
literature. For consistency, we re-determined Afinit and 
AiwD from the published Tcff and log g using our adopted 
WD and stellar evolutionary models. These clusters, the 
adopted parameters from each cluster, and references for 
the parameters and WD observations are given in Table 
[6l Points from the literature with small error bars typ- 
ically do not include the external fitting errors that we 
have considered in our data. 

As we are seeking the cleanest-possible IFMR in- 
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4 6 
Progenitor Mass (Mq 

Fig. 8. — The empirical initial-final mass relation. Points for 
M35 (filled circles) are shown along with cluster WDs from the lit- 
erature (triangles with gray error bars and lower limits) . Errors arc 
68% confidence limits and include only the propagated total fitting 
errors, not uncertainties in cluster ages. Curves are linear fi ts to the 
IFMR as determined by this work (solid line) and Fcrrario et all 
poos, dotted line), the semi-empirical IFMR of Weidemann (2000, 
dot-dashed curve), and the core mass at the end of thermal pulsing 
for Z = 0.008 and Z = 0.019 AG B stars (short- and long-dashed 
lines, respectively: IMarigo &: Girardiii2007i) . 



TABLE 6 
Adopted cluster parameters for 
previously-published white dwarfs 



Cluster 


Age (Myr) 


[Fe/H] 


References 


Pleiades 


125 ± 25 


0.00 


1,2,3,4,5 


NGC 2516 


158 ± 20 


-0.10 


6,7,8 


M34 


225 ± 25 


0.07 


9,10,11 


Sirius A/B 


238 ± 13 


0.00 


12,13 


NGC 2287 


243 ± 40 


0.00 


14 


NGC 3532 


300^25 


0.00 


6,14 


NGC 2099 


490 ± 70 


0.05 


14,15,16 


NGC 6633 


562^^9 


-0.10 


17,18 


Hyades 


625 ± 50 


0.13 


3,19,20 


Praesepe 


625 ± 50 


0.13 


3,21,22 


NGC 7789 


1400 ± 140 


-0.10 


23 


NGC 6819 


2500 ± 250 


-0.02 


23 


NGC 6791 


8500 ± 1000 


0.35 


24,25 



References. — (1 ) [Mazzci & Pigatto 198^, 
(2)IStauffer et al.lll998l, (3) Clavcr ct al. 2001, (4) 
[Ferrario ct al. 2005, (5) Dobbie ct al. 2006b, (6) 
IMev net ct al. 1993, (7) Koester & Reimers 1996, 
(8) teung ct al. 2002, (9) Jones & Prosser 199^, 
ISchuler et al., 2003. ( 11) .Rub in et al. 



.12) ILiebert et al.ri2005bl. (13) [Barstow e t al 
2003. f 14) [Dobbie et al.ll200a ( 15) [Marsh all et al, 
200E . (16)'Hartman ct al. 2008', (17) Jeffries et al, 
200a. (1 8) Williams & Boltc 20071, (19") 
Boeseaard & Friol 1990, (20) Ferryman e t al 



199a. (21) Dobbie ct al. 2004, (22) Dobbic c t al, 
2006al . (23) Kalirai ct al. 2008, (24) lOriglia et al, 
20061. (25) Kalirai et al. 2007 



dicative of the results of single-star evolution, a few 
WDs from the literature are excluded from our anal- 
ysis. WD 0437+138 is of spectral type DBA; the 
cooling models we employ are appropriate only for 
DA WDs. WD 0837-M99 is strongly magnetic, which 



could potentially i ndicate a more co mplicated evolution- 
ary history (e.g., iTout et all I2008D. T h e two binary 
WD candidates from IWilliams fc Boltd (|2007l) . NGC 
6633:LAWDS 4 and NGC 6633:LAWDS 7, either have 
uncertain parameters (if binaries, the observed spec- 
tra are a blend of the two components) or, if lone 
WDs, are not cluster members. Proper motion stud- 
ies have ide ntified the following WDs as field stars: WD 
0837-1-218 fC asewell et alll2008h . NGC 1039:LAWDS 9, 
NGC 1039:LAW DS 20, NGC 1039:LAWDS SI, and NGC 
7063:LAWDS 1 (iDobbie et al. 2008). 

From the figure, we see that the M35 WDs fall along 
the existing empirical IFMR, with higher cluster WD 
masses are seen to correspond to higher progenitor star 
masses, as is generally expected from stellar evolution. 
In fact, the M35 W Ds follow the linear relation of 
iFerrario et al.l ()2005D . perhaps with a shghtly steeper 
slope. Als o included in the fi gure are the semi-empirical 
IFMR of IWeidemaiml ()2000D and a theoret ic al pr edic- 
tion for the IFMR from ' Marigo fc Girardil (|2007f ) for 
Z = 0.008 and Z = 0.019 stellar evolutio nary models. 

A recent paper bv lSalaris et al.l (120081 ) highlights that 
the IFMR determined via our method is not completely 
self-consistent. The cluster age was adopted from litera- 
ture using isochrone-fitting methods employing different 
stellar evolutionary models; the WD evolutionary models 
do not use the chemical profiles output by the same stel- 
lar models used to determine the cluster age; and the at- 
mospheric models used to fit the spectra are not identical 
to those used in th e WD evolutionary models. However, 
as pointed out bv ISalaris et al.l (j2008l ). the overwhelm- 
ing source of systematic error in IFMR calculations is 
from the uncertainty in the cluster age (see §4.3.ip . The 
systematic errors due to the self-inconsistency described 
above is therefore assumed to be negligible. 

4.2.1. The shape of the IFMR 

The empirical IFMR can be fit marginally well by a 
linear function. We determine the slope and zero-point 
of this line using least-squares fitting with error bars in 
two dimensions (e.g., Pres s et al- 1992). We include the 
61 WDs in Figure [5] with masses Mwd > O.SM©; lower 
masses, which may indicate He cores and a different evo- 
lutionary history, were omitted. LAWDS 11 was also 
omitted from the fit due to its uncertain membership; if 
it is a cluster member and is the result of binary evo- 
lution, this exclusion would also be justified. Using the 
published random error bars for each WD's parameters, 
we obtain 

Mfinai = 0.339 ± 0.015 + (0.129 ± 0.004)Minit ; (4) 

this fit is plotted in Figure |5| The reduced value for 
this fit is 2.47. We note that, if LAWDS 5 and LAWDS 6 
are also excluded from the fit due to their poor spectral 
fits, the linear fit is unchanged within the formal error 
bars. 

Our linear fit is f ormally incons is tent w i th the slopes 
of the fits given in IFerrario et al.l (|2005f l. iKalirai et al.l 
(j2008l ) and ICatalan et all (j2008l ): the M35 WDs at the 

high-mass end of the relation prefer a steeper slope. How- 
ever, the goodness-of-fit statistic for all of these linear fits 
(including the new fits presented here) are very low. This 
could indicate that the error bars on the WD parameters 
are understated, that there is an additional scatter not 
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due to measurement error (such as intrinsic scatter in the 
relation), and/or that a linear functio n is not an appro - 
priate representation of the data fe.g.. lPress et al.iri992D . 
In any of these cases, the formal error bars on the slope 
and zero point of the linear fit quoted above are likely 
understated. Qualitatively, all of these linear fits appear 
to be a decent approximation to the data. 

We cannot exclude the possibility that the empiri- 
cal IFMR may not be strictly linear; emergent hints 
of a non-linear I FMR have be en discus sed recently in 
the literature fe.g.lKaiirai et al 112008; Salaris et al.ll2008l: 
[Catalan et aIll2008HDobbie et al.ll2008l) . The formal sig- 
nificance of these non-linear fits is limited by the larger 
error bars in massive WDs (driven by cluster age uncer- 
tainties and the steep change in the relationship between 
stellar mass and lifetimes at these younger ages) and the 
limited ar nount of data at the low Mjnit end of th e rela- 
tion (e.g.. [Catalan et al.ll2008t ISalaris et al.ll2008[ ). 

4.3. Systematic errors in the IFMR 
4.3.1. The age uncertainty in M35 

Although several systematic effects may affect the 
IFMR, by far the largest systematic in M35 is the un- 
certainty in the cluster age. Errors in the cluster age 
affect the IFMR in two distinct ways. First, an increase 
(decrease) in the assumed cluster age will result in all 
WDs from a given cluster having lower (higher) derived 
initial masses. Second, a change in the cluster age results 
in a larger shift for cluster WDs with higher progenitor 
masses, as the change in age represents a larger fraction 
of the main sequence lifetime for these stars. 

Both of these effects can be illustrated in the M35 em- 
pirical IFMR. Figure [5] shows the empirical IFMR for 
three different assumed ages of M35: 150 Myr, 175 Myr, 
and 200 Myr. The figure illustrates that, for older as- 
sumed ages, the cluster points shift systematically to 
lower Minit. The slope of the cluster IFMR also steepens 
significantly with increasing assumed cluster ages. 

The conclusions drawn from comparison of the M35 
IFMR to that of other clusters will differ depending on 
the assumed cluster age. If the cluster age is 150 Myr, 
the cluster IFMR agrees extraordinarily well with the 
empirical IFMR derived from other clusters. Yet if the 
cluster age is 200 Myr, one would conclude that the clus- 
ter WDs are systematically more massive than those in 
other young clusters, perhaps i ndicative of metallicity- 
dependent niass lo ss rates fe.g.. iMarigo fc Girardi|[2007t 
IKalirai et all[200l . 

These same systematics affect attempts to calculate 
the intrinsic scatter in the IFMR. For a younger M35, 
the scatter at a given Mjnit would be explained by the 
intrinsic observational errors. Yet for an older cluster, 
the intrinsic scatter would be significant. 

This systematic is most severe in young clusters; for old 
open clusters, the uncertainty in the cluster age is small 
relative to the progenitor star lifetimes. For example, 
the 1 Gyr (- 15%) uncertainty in the age of NGC 6791 
results in just a 0.05Afo change in Mjnit for its WDs, 
compared with the > IM© uncertainty in Mj^it for the 
M35 WDs for the cluster age uncertainty of just 25 Myr 
(also - 15%). 

4.4. Limits on M^r, 



Observational lower limits on the maximum mass of 
WD progenitors, Mup, can be determined from IFMR. 
To zcroth order, this can be accomplished by simply iden- 
tifying the WD with the highest Afinit in the empirical 
IFMR. However, this approach is complicated by numer- 
ous uncertainties. First, the sizes of the error bars on in- 
dividual WDs are large enough that it is not possible to 
identify the single WD with the highest Minit. Second, 
the errors on Minit of individual WDs are asymmetric 
and non-analytic, complicating any effort to simply cal- 
culate limits on A/up- Third, the WDs with the highest 
Minit come from numerous star clusters, each of which 
will have its own unique systematic errors in Minit due to 
the cluster age uncertainties. We therefore elect to use a 
Monte Carlo simulation to determine the lower observa- 
tional limits on A/up- 

The simulation creates 25000 realizations of the ob- 
served WD population in the open clusters M35, 
the Pleiades (inclu ding CD 50 and PC 0136-^251, 
IDobbie et al] l2006H). NGC 2516, NGC 3532, NGC 2287, 
and Sirius B; WDs from older open clusters are excluded 
because they either have such low Minit or such large 
error bars that they do not contribute to the limits on 
Afup. For each realization, each open cluster is a as- 
signed an age randomly drawn from the values quoted in 
Table [51 as the published errors in cluster ages tend to 
indicate acceptable ranges rather than a normal distri- 
bution about the best value, cluster ages are drawn ran- 
domly from a uniform distribution. Additionally, in each 
realization a T^ff and logg are assigned to each known 
WD, with values drawn from a normal distribution with 
a mean of the observed Toff and log 5 for the WD, and a 
standard deviation equal to the stated random errors in 
each quantity. The assigned Tes and log 17 are then used 
to determine Minit for each WD; if the cooling age of the 
simulated WD is older than the cluster age, that WD is 
ignored for that particular realization. For each realiza- 
tion, the highest overall Minit is identified as the lower 
limit on Mup for that realization. 

We perform six runs on this Monte Carlo simulation 
with varying parameters. Three different WD popula- 
tions are tested: 

• Only the WDs in NGC 2168; this limits the num- 
ber of systematics involved in the simulation, but 
provides the loosest limits on Afup, 

• All DA WDs in the five open clusters mentioned 
in the previous paragraph as well as Sirius B; this 
data set represents all available useful cluster data, 
and 

• The previous sample minus three WDs of uncer- 
tain cluster membership: NGC 2168 LAWDS 11 
(see SH), and GD 50 and PC 0136-1-251, as their 
common space motion with the Pleiades does not 
prove they were born simultaneously with the clus- 
ter. 

Each of these three WD populations was paired with two 
variations on WD core composition. In the first instance, 
all WDs are assumed to have carbon-oxygen cores; in the 
second, any WDs with Mwd > 1.O5M0 are assumed to 
have oxygen- neon cores (see H4.6.ip . 

Results of these calculations are presented in Table [71 
Based on the WDs in M35 alone, the 95% confidence 
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Fig. 9. — The effects of assumed cluster age on the empirical initial-final mass relation. The M35 points shift systematically to lower 
Minit and their slope steepens as the assumed cluster age is increased from (left) 150 Myr to (center) 175 Myr and (right) 200 Myr. Symbols 
are as in Figure |8] 



TABLE 7 

Lower Limits on Mi 



WD Sample Carbon/Oxygen Oxygen/Neon 

50% 90% 95% 99% 50% 90% 95% 99% 
(Mq) (Mq) (Mq) (Mq) (Mq) (Mq) (Mo) 



M35 alone 6.57 5.40 5.19 4.89 6.48 5.34 5.14 4.86 

All Clusters 8.86 7.41 7.08 6.51 8.39 6.95 6.70 6.34 

All Clusters, "Cleaned" 7.97 6.58 6.25 5.77 7.80 6.54 6.30 5.91 



Note. — Masses are lower limits on Mup from Monte Carlo simulations for the given 
confidence interval and given core composition for 1^ 1.05Mq (lower mass WDs are 

assumed to have C/O cores). The "cleaned" sample excludes LAWDS 11, GD 50, and PC 
0136 + 251 



lower limit on Mup is 5.1Mq — 5.2Mq, depending on 
assumed core composition. By combining all available 
cluster WD data, the 95% confidence lower limit on A/up 
increases to 6.3Mq — 7.1Mq, depending on the member- 
ship of LAWDS 11, GD 50, and PG 0136-)-251 and on 
the assumed core composition. 

We also note that the combination of all cluster data 
is required to get the tightest limits on Mup- For exam- 
ple, the Pleiad WD, LB 1497, has the highest plotted 
Afinit in Fig. [H In our Monte Carlo calculations exclud- 
ing LAWDS 11, GD 50, and PG 0136-1-251, LB 1497 
only provides the highest Afinit in 38% of the realiza- 
tions. Messier 35 WDs provide the highest Afinit in 42% 
of the realizations, and NGC 2516 in 19% of the real- 
izations. Excluding any one of these three clusters from 
the Monte Carlo calculation lowers the limit on Afup by 
~ O.5M0 - lAf0. 

Combining these lower limits on Afup with emerging 
upper limits from su pernov a prog enitor searches (Afu p < 
9.5Af0, e.g.. iHendrv et al.ll2006riSmartt et al]l2008f ). it 
may be reasonable to assert that 6Mq < Mup < 9.5M0 
with 95% confidence for stars near solar-metallicity. 
However, a more careful analysis of both the WD and 
supernova data is needed before these constraints can be 
claimed with confidence. 

4.4.1. Systematic errors in the limits on Mup 

The primary sources of systematic error in these simu- 
lations are from the input star cluster ages and the input 
stellar evolutionary models. The input ages are all de- 
termined from evolutionary models including moderate 
convective overshoot, and so are reasonably consistently 



determined despite coming from a multitude of published 
sources and methods. However, there is a well-known de- 
generacy between an observed cluster's age, metallicity, 
and distance; a significant change in the age of any of the 
clusters used in this simulation would change the results. 
And, as stellar models continue to evolve, the derived 
ages for the entire cluster sample may change systemat- 
ically. 

These same evolutionary models are used to con- 
vert the WD progenitor lifetime into Minit- There- 
fore, any significant changes in the output lifetimes for 
intermediate-mass stars would also significantly impact 
these results. However, moderate levels of convective 
overshoot have become generally accepted in the com- 
munity and appear to be borne out by o bservations (e.g., 
iClaretl [20071 : iSandberg Lacv et al.l[2008l ) . 

How large a difference could future changes to evo- 
lutionary models make? According to the Padova 
isochrones, a solar-metallicity model including moderate 
convective overshoot with a lifetime of 100 Myr has a 
mass of 5.42M0, and a mass of 5.18Af0 with no con- 
vective overshoot. Therefore, it seems unlikely that our 
derived limits will change by < 0.3 A/© due to a choice of 
evolutionary models alone. Larger changes are possible, 
though, depending on how newer evolutionary models 
affect the input star cluster ages. 

Another possible source of error is due to our combi- 
nation of data from star clusters of differing metallicity; 
namely, NGC 2168 is significantly subsolar metallicity, 
while the remaining clusters and Sirius are close to so- 
lar metallicity. Numerous models of super-AGB stars 
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exhib it a dependence of M^p on metallicity (e.g., ISiessI 
l2007l and references therein). Though these models dif- 
fer significantly on the precise value of Mup, they pre- 
dict that Mup for stars with the metallicity of NGC 2168 
should be systematically ~ 0.2Mq to 0.3Mq lower than 
for solar-metallicity stars. 

We note that incompleteness in the WD sample cannot 
reduce our lower limits on Mup; additional WDs can only 
raise the lower limit. 

4.5. The cluster white dwarf distance modulus 

The distance modulus and reddening we have adopted 
for M35 [{m-M)v = 10.3±0.1, E{B-V) = 0.22±0.03] 
are based solely on previously-published main-sequence 
fitting derivations. As the spectral fits to each WD, when 
combined with the WD photometry, provide a measure 
on the distance modulus and reddening for each WD, we 
can use the ensemble averages of these individual mea- 
surements as an independent determination of the cluster 
distance and reddening. Again, due to the uncertainty 
in its spectral parameters, we do not include LAWDS 22 
in these calculations. 

The weighted means of the individual distance mod- 
uli give a mean WD distance modulus of (m — M)v = 
10.45 ± 0.08, with a dispersion of (J(^rn-M) — 0.28 mag. 
The mean color excesses are E{B — V) — 0.185 ± 0.010, 
aE(B-v) = 0.055 and E{U -V) = 0.329 ± 0.014, 
'^E(u-v) = 0.091. The mean cluster WD distance mod- 
ulus and E{B — V) are therefore fully consistent with the 
values determined from main-sequence fitting; the ratio 
E{U-V)IE{B-V) = 1.78 ± 0.12 is consisten t with the 
Milky Way reddening law of lRieke fc Lebofskvl ()1985[ ) for 
Rv = 3.1. 

The measured scatter in the measured color excesses is 
larger than the uncertainties in the individual measure- 
ments (see Table [5]) . In addition, there is a significant 
correlation between the individual WD distance moduli 
and the E{B~V) color excess. These two effects could 
be explained by differential interstellar extinction across 
the cluster. If we "correct" each distance modulus to 
the mean E{B — V) color excess, the dispersion in the 
distance moduli drops to cr(^,n-M) — 0-19 mag. 

However, there is no obvious spatial correlation be- 
tween WDs with similar measured E{B — V). In addi- 
tion, errors in the measured Tcs will mimic extinction 
effects. If the measured Teg is higher than the actual 
value, then the derived WD absolute magnitude will be 
too bright, resulting in too large of a measured distance 
modulus. At the same time, the derived model color will 
be bluer than the actual WD, resulting in an artificially 
large measured E{B — V). The magnitude of this effect is 
essentially identical to the stated errors in the E{B — V) 
and E{U — V) measurements in Table [5l Therefore, the 
errors in Teg can explain most, but not all, of the ob- 
served scatter in the E{B — V) and E{U — V) values. 

4.6. Notes on interesting objects 

4.6.1. Potential oxygen-neon core WDs 

Four M35 WDs, LAWDS 2, LAWDS 14, LAWDS 27, 
and LAWDS 30, have masses within la of 1.05Mo, the 
lower mass of ON e core WDs that ma y be produced by 
super- AGB stars ()Althaus et al.ll2007D . For this reason, 
we have also calculated Mwd and Mjnit for these stars 



using the ONe WD evolutionary models of lAlthaus et al.l 
(|2007( ). The ONe models reduce Mjnit by ^ 0.2Mq for 
all four WDs. The decrease in Mi„it is due primarily 
to the lower heat capacity of ONe cores, which allow the 
ONe-core WDs to cool more rapidly than C/O-core WDs 
of the same mass. 

From our data, we have no means of knowing if these 
two WDs have carbon-oxygen or oxygen-neon cores. 
However, stellar evolutionary models of the ~ 5Mq pro- 
genitors of these WDs suggest that they should have C/0 
cores. 

4.6.2. The DB white dwarf LAWDS 4 

Claims have been made that a lack of DB WDs 
exists in younger open clusters (|Kalirai et al.l l2005al 
[2008h . Several non-DA WDs are known to exist in 
the field of younger open clusters, including LP 475- 
252 (spectral type DBA in the Hyade s), NGC 2168: 
LAW DS 28 (the hot DQ in this cluster; IWilliams et"an 
[2006h. NGC 1039:LAWDS 26, and NGC 6633 :LAWDS 
16 f Willia ms fc Boltd [20071: iRubin et"an[200l . though 
only LP 475—252 has been confirmed as a cluster mem- 
ber via proper motion. We have identified one DB WD 
in the field of M35, LAWDS 4. If this is a cluster mem- 
ber, it would make M35 the youngest open star cluster 
with a DB WD member (noting again that M35 contains 
the non-DA LAWDS 28). 

Is LAWDS 4 a cluster member? At present, the 
low signal-to-noise of the spectrum precludes a robust 
spectral fit. However, we can estimate its T^ff. The 
weakness of the helium absorption lines, the spectral 
slope, and the colors are all consistent with a fairly cool 
(Toff w 15000- 17000K) DB WD. Based on coohng mod- 
els provided by P. Bergeron, this DB could only be this 
cool and still be younger than M35 if it has a relatively 
low mass (M < O.eM©). 

If we assume a WD mass of 0.4Afo to O.6M0 and Toff — 
17000K, the apparent distance modulus of LAWDS 4 is 
between 10.24 and 10.77, consistent with the cluster dis- 
tance modulus. In other words, this DB is photometri- 
cally consistent with cluster membership, under the as- 
sumption that its mass is M < O.6M0. 

If LAWDS 4 is a cluster member, its progenitor mass 
would have to be high (> 5Mq), yet the WD mass 
(< O.6M0) lies well below the empirical initial-final mass 
relation, possible under a binary formation scenario for 
the WD. It is also possible, and perhaps probable, that 
this star is a field WD located near the cluster. Accu- 
rate, precise proper motion measurements for this WD 
will likely be necessary to clarify its cluster membership, 
and a higher signal-to-noise spectrum is needed to deter- 
mine its Toff and log 5. Clearly, this star warrants further 
study. 

5. CONCLUSIONS 

In this paper, we have presented observations and anal- 
ysis of the white dwarf population of the open star cluster 
M35. Our conclusions arc the following. 

• Spectroscopy of 14 white dwarf candidates identi- 
fies 12 DA, 1 DB, and 1 hot DQ WDs. Temper- 
atures, surface gravities, masses, and cooling ages 
are derived for each of the DA WDs. 
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• All twelve DAs are potentially cluster member 
WDs, with distance moduli consistent with that 
of the cluster and cooling ages less than the clus- 
ter age. The hot DQ and DB are also consistent 
with being cluster members. Further data, such 
as proper motion measurements, are necessary to 
confirm the cluster membership of each star. 

• The empirical initial-final mass relation from the 
M35 white dwarf population is consistent with the 
roughly linear relation derived from other open star 
clusters. 

• The dominant systematic uncertainty in the em- 
pirical initial-final mass relation of M35 is the un- 
certainty in the star cluster age. In the absence 
of tighter age constraints, we cannot draw robust 
conclusions on the intrinsic scatter and metallicity- 
dependence of the initial-final mass relation. 

• Based on M35 WDs alone, the lower limit on the 
maximum mass of WD progenitor stars (M^p) is 
~ 5.IM0 (95% confidence). Inclusion of WDs from 
other young open clusters raises this lower limit 
to ~ 6.3M0 — 7.IM0, depending on the member- 
ship of certain massive white dwarfs and the core 
composition of the most massive WDs. Combined 
with upper limits on Mup from supernova surveys, 
QMq < Mup < 9.5Mq for solar-metallicity stars. 

• Based on the cluster WDs alone, we derive a dis- 
tance modulus to M35 oi{m-M)v = 10.45 ± 0.08 

and reddening of E{B-V) = 0.185 ± 0.010, both 
in agreement with published values derived from 
main sequence fitting. 

• Four of the DAs have masses that are sufficiently 
massive that they may possess oxygen-neon cores. 
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